Inhomogeneous metallic phase upon disordering a two dimensional Mott insulator 
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We find that isoelectronic disorder destroys the spectral gap in a Mott-Hubbard insulator in 2D 
leading, most unexpectedly, to a new metallic phase. This phase is spatially inhomogeneous with 
metallic behavior coexisting with antiferromagnetic long range order. Even though the Mott gap 
in the pure system is much larger than antiferromagnetic exchange, the spectral gap is destroyed 
locally in regions where the disorder potential is high enough to overcome the inter-electron repulsion 
thereby generating puddles where charge fluctuations are enhanced. With increasing disorder, these 
puddles expand and concomitantly the states at the Fermi energy get extended leading to a metallic 
phase. We discuss the implications of our results for experiments. 



Disorder can lead to new importantphenomena: such 
as localization in disordered media [ij, |2|, metal-insulator 
transitions 0, precisely quantized Hall plateaus in the 
quantum Hall effect, and zero resistance from pinned 
vortices in a type II superconductor. In systems as 
diverse as superconducting films0], high T c cuprates[^ 
and manganitesQ disorder produces nanoscale inhomo- 
geneities which can have a profound influence on the 
properties of these systems. 

We propose the existence of a new metallic phase at 
T = which is generated when potential disorder is in- 
troduced into a Mott insulator. This metallic phase is 
highly inhomogeneous and coexists with long range anti- 
ferromagnetic order. It is sandwiched between two qual- 
itatively different insulators: at low disorder, a Mott in- 
sulator with a spectral gap and antiferromagnetic (AFM) 
long range order, and, at high disorder, a paramagnetic 
Anderson insulator with gapless excitations. Our result 
is surprising since one would expect that the primary role 
of disorder is to localize wave functions, which is opposite 
to the conditions for forming extended states necessary 
for metallic behavior. 

We find that the route to the formation of a metal at 
intermediate disorder is the following: Given a random 
potential profile, there are two types of regions: weakly 
disordered (WD) with relatively small potential fluctua- 
tions and strongly disordered (SD) with larger fluctua- 
tions. AFM order exists in the WD regions that consist 
of singly occupied sites with a large local Mott gap. On 
the other hand, in the SD regions the local density devi- 
ates from one per site which disrupts the antiferromag- 
netic ordering and also leads to the collapse of the gap 
locally. There is thus a 'phase separation' between the 
insulating AFM WD regions and the metallic SD regions 
with suppressed AFM. The low lying excitations of the 
system remain trapped in the SD regions. With increas- 
ing disorder the SD regions grow and along with that, 
the eigenfunctions get extended throughout the system 
and the system becomes metallic at intermediate disor- 
der. At high disorder the system becomes a paramagnetic 
Anderson insulator. We substantiate the above picture 



with calculations of the magnetic correlations, the local 
density of states and the optical conductivity. 

We propose that experiments on cuprates in their par- 
ent Mott insulating phase, could introduce disorder by 
isoelectronic substitution, or by bombardment, without 
changing the carrier doping. Such experiments should be 
able to see this unusual metallic phase at intermediate 
disorder. Recently, Vajk et. al(2| have measured the spin 
correlations using neutron scattering in LC12CUO4 with 
Zn substitution on the Cu sites. It will be necessary 
to complement the magnetic information from neutron 
scattering with transport and spectroscopy to look for 
signatures of an inhomogeneous metallic phase. 
Model and Calculation: We model the 2D disordered 
Mott insulator by a repulsive Hubbard model at half fill- 
ing (one electron per site) with potential disorder: 



H = K + n iT n tl + /X V i ~ V) n i<y 



(1) 



K = -iE<y>, CT ( c L c jV + h.c) is the kinetic energy, c] a 
(cicr) the creation (destruction) operator for an electron 
with spin a on a site Yi of a square lattice of N sites, 
t the near-neighbor hopping, U is the on-site repulsion 
between electrons, — c\ a ac an d n the chemical po- 
tential. The random potential Vi is chosen independently 
at each rj and is uniformly distributed in the interval 
[—V, V]; V thus controls the strength of the disorder. All 
energies are measured in units of t. This is a minimal 
model containing the interplay of electronic correlations 
and localization: for zero disorder V = Q it describes a 
Mott insulator at half filling, and for U = it reduces 
to the (non-interacting) Anderson localization problem. 
We report results for U = 4t and lattices up to 50 x 50. 

The dynamical mean field theory has been widely 
used to study Mott insulators and includes the quantum 
fluctuations associated with an effectively single site ap- 
proximation of the system. The inhomogeneous Hartree 
Fock approximation that we implement here is a comple- 
mentary method which does not contain the dynamical 
fluctuations but, when solved self consistently as we do 
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FIG. 1: (a) Left panels show the local staggered magneti- 
zation m) {i) = (S z (i))( — iy* +ty . The regions in blue have 
AFM order; the defective sites with reduced AFM order are 
shown in yellow and red. The data is for the disordered Hub- 
bard model at half filling with U = At on a 28 x 28 lattice for 
disorder strengths V = 2t; V — 3t; V = 5t for one realization 
of disorder, (b) Probability distribution P(m t ) of m'(i) for 
different values of V. For V = 1, P(m t ) has a peak near 0.4. 
With increasing V, P(mJ) gets broader and develops weight 
near indicating the growth of paramagnetic regions, (c) 
Probability distribution P(n) of the site-occupancy (m) for 
different values of V. For V = 1, P{n) has a peak near an 
average single site occupancy (m) ft! 1. With increasing V, 
P(n) gets broader and develops weight for doubly occupied 
and unoccupied sites. 



here, captures the inhomogeneity in the spatial struc- 
ture of the electronic density, a feature which is at the 
heart of problems with disorder. The only known way 
to include both the spatial and temporal fluctuations is 
quantum Monte Carlo (QMC) but that method too has 
limitations arising from finite size effects, inability to ex- 
tract dynamical information without the use of compli- 
cated maximum entropy methods and the sign problem 
for fermions that limits the calculations to finite tem- 
peratures. We begin by treating the spatial fluctuations 
of the local densities riia = (cj^Cjo-) and local magnetic 
fields hf = —U (cj|Cij) and = —U (c|, Cjf) using a site 
dependent mean field approximation. As a consequence, 
we get an effective Hamiltonian which is quadratic in the 
fcrmion operators and can be diagonalized for a system 
with periodic boundary conditions. Starting with an ini- 
tial guess for hf's and ra^'sjl we numerically solve for 
the eigenvalues e„ and eigenvectors ^w(rj) of a 2N x 2N 
matrix. The local fields are then determined in terms of 
the eigenfunctions and the eigenvalue problem with these 
new local fields as input is then iterated until self consis- 
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FIG. 2: (a) Density of states averaged over 10 realizations at 
half filling and U = At for V/t =1, 2, 3. For V = t there is 
a Mott gap in the spectrum at the chemical potential E — 0; 
and for V — 2t the gap has closed, (b) Single particle energy 
gap as a function of V showing the collapse of the gap at 
Vci ~ 2.2t (scale on LHS). The decrease of the AFM staggered 
order parameter m) (scale on RHS), as a function of disorder 
strength V and its vanishing beyond V C 2 ~ it > V c i- Note 
that at V — 2t the system is still strongly AFM even though 
the Mott gap has vanished. 



tency is achieved at each site. 

Antiferromagnetic Order and Spectral Gap: Fig. 1 
shows that the local staggered magnetization rtv(i) = 
(— l) lx+i v (S z (i)} is largely antiferromagnetic for V = 2t 
with a few local defective regions. With increasing dis- 
order the defective regions having reduced antiferromag- 
netic order grow in size. 

We next look at the behavior of the two defining char- 
acteristics of a Mott insulator in Fig. 2: the spectral 
gap obtained from the lowest eigenvalue of the effective 
Hamiltonian, and the AFM order parameter m'(i) ob- 
tained from the spatial behavior of the spin-spin corre- 
lation function at large distances, as a function of dis- 
order. The first surprise is that even though the energy 
scale for charge fluctuations is U 3> J ~ t 2 /U, the scale 
for antiferromagnetic exchange, the spectral gap vanishes 
at Vci ~ 2.2i which is lower than the critical disorder 
V C 2 ~ At where AFM long range order (AFLRO) van- 
ishes. 

The occurrence of two critical disorder strengths de- 
fines three distinct regions: In region I defined by < 
V < V c i the system is a Mott insulator with a finite gap 
and AFLRO. The intermediate region II, V c \ < V < V c2 
is extremely unusual with AFLRO but no gap and is dis- 
cussed in greater detail below. In region III, V > V C 2 the 
system is simply an Anderson or localized insulator with 
gapless excitations and no magnetic order. 
Nature of eigenfunctions: We argue below that in re- 
gion II the system is a disordered AFM metal. How 
does this metal come about? The spatial extent of 
the eigenstates near the Fermi energy in Fig. 3(a) dis- 
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FIG. 3: (a) Single particle eigenstates at the Fermi energy 
for disorder strengths V = 2t, 3t, 5t. Note that the eigen- 
states are localized for low V = 2t and high disorder V = 5t 
but surprisingly are more extended at intermediate disorder 
of V = 3t. (b) Inverse participation ratio (IPR) or the local- 
ization length for the wave function at the Fermi energy 
as a function of disorder V. In both regions I and III IPR 
increases with disorder or £; oc decreases with disorder, as ex- 
pected. But in the intermediate region II, the wave function 
in fact gets less localized with increasing disorder. 



tinctly shows that they are localized at low V = 2t 
and high disorder V — 5i, but quite surprisingly at 
intermediate disorder V = 3i the states get more ex- 
tended. The localization length £i oc (aO for the (normal- 
ized) state Tp a (ri) is related to the inverse participation 
ratio IPR = J2 r \ip a (ri)\ 4 oc £ioc( a ) O ne would ex- 
pect that with increasing disorder, £i oc would decrease or 
correspondingly IPR would increase. Instead we see in 
Fig. 3(b) a very definite decrease in IPR for V ~ U/2 sig- 
naling that even though disorder is increasing the states 
are getting more extended. Such anomalous behavior 
continues till about V ~ 3J7/4 and then once again re- 
verts to the usual behavior where IPR increases with V . 
Based on numerical calculations on finite size systems it 
is difficult to prove if the states in region II are truly 
extended, but what is abundantly clear is that the IPR 
does show a very definite non-monotonic behavior signal- 
ing a marked change in the nature of the states, at least 
on mesoscopic scales. 

Previous QMC simulations on disordered bosonic[ll| 
and fermionic|l2| have shown that potential disorder that 
breaks particle-hole symmetry produces superfluid stiff- 
ness or conductivity which can increase with disorder. At 
a Hartree Fock level such an effect can be understood as 
a screening of the random potential by the interactions 




FIG. 4: Frequency dependence of the conductivity for V — It 
in region I showing a gap uuo in the joint density of states; for 
V = 3t in region II showing a linear dependence indicative 
of metallic behavior; and for V = 5t in region III showing 
ui 3 dependence indicative of Anderson localization. The inset 
shows the decrease of the gap u>o with increasing disorder V 
and its vanishing around V c i = 2.2t. 



between the particles|l5j. Our calculations go beyond 
simple screening arguments and have revealed the phys- 
ical origin of the anomalous behavior seen in Fig. 3 at a 
microscopic level. 

Frequency- dependent conductivity: The nature of the fre- 
quency dependent conductivity Rea(co) = ImA(cj)/w 
gives insight into the conducting properties of the phases. 
A(u/) is the Fourier transform of A(r) = (j(t)j(O)), 
the (disorder averaged) current-current correlation func- 
tion. As shown in Fig. 4 the low frequency behavior of 
uj<j(uj) has a finite gap in the Mott region I for V = t; 
shows a linear uj dependence in the metallic region II for 
V = 3t which implies that there is a finite dc conductivity 
cr(u> — > 0); and shows a ui 3 dependence in the Anderson 
insulating regime III for V = 5f, which implies that there 
are gapless excitations but nevertheless a — > 0. 
Origin of the Metallic Phase: Now that we have seen 
evidence for a metallic phase at intermediate disorder, 
the question arises as to how such a phase comes about. 

In the clean system, the repulsive interactions prevent 
double occupancy. But once there is potential disorder, 
electrons can hop into sites that have deep potential wells 
and thereby lower their energy. This leads to a broadened 
distribution of local occupancy and a growth of paramag- 
netic sites, as illustrated in Fig. 1. In the classical limit, 
sites with a large positive potential Vi > U/2 are unoc- 
cupied; sites with —U/2 < Vi < U/2 are singly occupied; 
and sites with large negative Vi < —U/2 are doubly occu- 
pied. Only sites with one electron can contribute to AFM 
order if they form a connected cluster. The critical prob- 
ability p c for classical percolation of vacancies on a 2D 
square lattice is p c — 0.41 which should equal the density 
of doubly and unoccupied sites, so that p c = 1 — UJ (2V) 
which implies that (V/U) c = 1/2(1 -p c ) = 0.85. This 
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FIG. 5: (a) Histogram P(mr) of the local staggered mag- 
netization (black). We will define those sites with staggered 
magnetization 0.38 and more to be AFM sites (region shown 
in blue); and sites with (S z (i)) in the range (-0.1,0.1) to be 
paramagnetic PM sites (shown in red), (b) Histogram P(V) 
of the local disorder potential for all sites (black); histogram 
of AFM sites (blue) which coincide with the less disordered 
WD regions; and histogram of PM sites (red) in the SD re- 
gions, (c) Histogram P(n) of the local density for all sites 
(black); for AFM sites (blue); and for PM sites (red). 



indicates that an estimate of V C 2 ~ 3At for U = At. 
We have compared our results for as a function of V 
with the percolation of vacancies in a quantum Heisen- 
berg model and find very good agreement 0, 0] . 

Fig. 5, which shows the correlation between the 
strength of disorder and the nature of the regions, is 
very revealing. It clearly shows a reorganization or a 
phase separation of the system into WD regions that are 
AFM and have a local density of one electron per site, 
and the SD regions which are PM with a bimodal den- 
sity clustering towards zero and two per site. Our results 
for the existence of metallic behavior in 2D is even more 
surprising when we recall the scaling theory of localiza- 
tion for non- interacting electrons Q which claims that in 
dimensions less than and equal to two all the single par- 
ticle states are localized by arbitrarily small amounts of 
disorder. Thus the underlying picture that emerges for 
the origin of the unusual metal at intermediate disor- 
der in a Mott insulator at half filling is a percolation of 
the defected SD regions with reduced AFM. As seen in 
Fig. 3(a) (top panel) the lowest excitations live in these 
defected SD regions and with increasing disorder the de- 
fective regions expand and along with that the excitations 
get more delocalized as they are able to travel on the ex- 



panding SD clusters, see Fig. 3(a) (middle panel). This 
is a rather remarkable situation that increasing disorder 
in a correlated system is able to induce metallic behav- 
ior. With further increase of disorder the states on the 
percolating cluster get localized and the system becomes 
an Anderson insulator, see Fig. 3(a) (bottom panel). 
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